Multiple Changes of Order of the Vortex Melting Transition in Bi 2 Sr 2 CaCu208 with 

Dilute Columnar Defects 
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A low concentration of columnar defects (CDs) is reported to transform a first-order vortex lattice 
melting line in Bi2Sr2CaCu2 08 crystals into alternating segments of first- and second-order transi- 
tions separated by two critical points. As the density of CDs is increased, the critical points shift 
apart and the range of the intermediate second-order transition expands. A third, low temperature 
critical point was also observed in one sample. The measurement of equilibrium magnetization and 
the mapping of the melting line down to 27K was made possible by employment of the shaking 
technique. 
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PACS numbers: 74.25.Qt,74.25Dw, 74.72Hs, 64.70dg 
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The interplay between different energy scales of a phys- 
ical system can induce phase transitions and yield a rich 
phase diagram. In particular, in a clean system a melting 
transition results from the competition between elastic- 
ity and thermal fluctuations and is usually of a first-order 
nature. The presence of a third competing scale due to 
disorder can suppress it to second order (SO) and lead 
to an even richer picture, the understanding of which is 
still quite limited [H, ln order to gain insight into 
the various competing mechanisms, we introduce an ad- 
ditional control parameter in the form of a variable weak 
correlated disorder. We report an exceptionally intricate 
case, where the initially first-order (FO) melting transi- 
tion of the vortex matter in a high temperature supercon- 
ductor exhibits a rare three-section FO-SO-FO behavior 
and even seems to undergo a four-section FO-SO-FO-SO 
sequence. Our findings reveal that a SO transition is 
nucleated in a small segment along the FO transition 
line, bound by two critical points. This SO segment 
grows when the control parameter is increased thereby 
providing important insight into the general mechanism 
of transformation of a FO transition into a SO transition. 

Bi 2 Sr 2 CaCu208 (BSCCO) high-T c superconductor has 
an especially rich vortex matter phase diagram which 
results from the interplay between the elasticity, ther- 
mal fluctuations, and point disorder. It consists of a FO 
melting line Q which separates the low-field Bragg glass 
from the high-field disordered phases 0, H, S] • At high 
temperatures the FO melting transition is mainly ther- 
mally driven, whereas at lower temperatures it gradu- 
ally changes its nature into a disorder-driven FO inverse 
melting transition 0, IE El ■ Indications of an additional 
almost vertical SO glass transition line at intermediate 
temperatures were also found [IE ■ 

This phase diagram changes entirely with introduction 
of a fourth energy scale due to dense columnar defects 
(CDs) > B m (T), where B$ = n cd (j) , n cd is the 
density of CDs, <po is the flux quantum, and B m (T) is 
the temperature-dependent melting field. Compared to 



the correlated pinning energy, the elastic energy becomes 
negligible. As a result, the ordered solid phase is replaced 
by a disordered Bose glass (BoG) phase in which vortices 
are localized on CDs and melt through a SO BoG tran- 
sition [B, Ell- 
in this paper we focus on the dilute CD limit, where 
vortices outnumber CDs at most relevant fields. Re- 
markably, in this limit the energy due to correlated pin- 
ning is comparable to the other three energy scales. 
Consequently, the four relevant energy scales mold to- 
gether a particularly complex B — T phase diagram 
whose nature is not well-understood. Most experimen- 
tal [H, EE El, El [pL theoretical [U El, El and nu- 
merical studies [2fj, |2lj, [22 1 agree that at high temper- 
atures the transition remains FO. At intermediate tem- 
peratures where the melting occurs at intermediate fields 
B m > Bcf,, the pristine Bragg glass phase is replaced 
by a porous solid with the majority of CDs occupied 
by strongly pinned vortices that form an amorphous 
matrix. The remaining vortices form ordered crystal- 
lites that are embeded in the pores of the rigid matrix 
H EE EE EE EE EE HE ED, El ■ These vortex nanocrys- 



tals melt into nanodroplcts apparently through a second 
order transition. 

Experimentally, even a low density of CDs enhances 
hysteretic effects throughout most of the solid phase in 
bscco [IE El- This is a major obstacle in map- 
ping the true thermodynamic phase transitions of vor- 
tex matter. Lacking access to the thermodynamic be- 
havior, past ex per imental studies have focused on dy- 
namics 112, 13, 1 4 , Il6l, |23|1 . It was concluded that the 



high temperature FO transition line terminates at a crit- 
ical point (CP) whose location depends on CD density 




151 116H . consistent with theory and simulations 
IE EJ, ED, HI • At temperatures below the CP, 
however, these out-of-equilibrium studies could not eval- 
uate the exact nature of the transition. Thus, a SO tran- 
sition was premised to exist at all temperatures below 
the CP [IE El EE OS EJ ED, HI], though other theoret- 
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ical scenarios were also considered [17(. This incomplete 
understanding of the dilute pinning limit led naturally 
to the assumption that as CD density is increased, the 
SO nature of melting spreads from low temperatures to 
higher ones with a shift of a single CP. Our findings, 
however, suggest an essentially different process. 

Within the present study we mapped for the first time 
the equilibrium phase diagram of vortex matter with di- 
lute CDs by performing local magnetization measure- 
ments during vortex shaking. This method 

0, EH EI El 

utilizes an in-plane ac shaking field which agitates the 
vortices and assists them to assume their equilibrium po- 
sitions. Tilting the magnetic field away from the CDs was 
shown to alter the equilibrium magnetization in YBCO 
crystals [26|. In contrast, in BSCCO moderate in-plane 
field has essentially no effect on the equilibrium prop- 
erties even in pre sence of CDs due to the very high 
anisotropy (lR \2l\. As a result, the ac in-plane field 
enhances vortex relaxation without altering the thermo- 
dynamic transitions We measured three BSCCO 
samples irradiated at GANIL with matching fields of 
= 5G, 10G, 20G and critical temperature T c w 92K, 
and an additional sample with B$ = 20G and T c ps 91K. 
Half of each sample was masked during irradiation to al- 
low a direct comparison between pristine and irradiated 
behavior. The local magnetic induction of the sample 
was measured by an array of Hall sensors 10x10 /im 2 
each, while sweeping the external magnetic field H at 
a constant temperature. Simultaneous measurements of 
the transition on the pristine halves of the samples show 
that it remained FO throughout the temperature range. 
Shaking fields up to 350Oe at 10Hz were used. 

Figure [T] presents the measured local magnetization, 
B - H, at (a) 90K, (b) 84K and (c) 42K. A FO transi- 
tion appears as a sharp step in B(H) and a SO transition 
is manifested by a break in its slope. To better resolve 
the nature and location of the phase transition we differ- 
entiated the measured induction B with respect to the 
applied field H. Figure [2] shows the derivatives dB/clH 
of representative (a) high temperature, (b) intermediate 
temperature and (c) low temperature measurements. A 
(5-likc peak in the derivative dB /dH indicates a FO tran- 
sition, whereas a discontinuity signifies a SO transition. 

At high temperatures we find that the transition re- 
mains FO in the presence of CDs, similarly to that in the 
pristine areas, as previously reported [13l . 
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FIG. 1: (color online) The measured local magnetization, 
B — H , of the = 10G sample in the presence of shak- 
ing, measured upon ascending and descending external field 
at constant temperatures of (a) 90K, (b) 84K and (c) 42K. A 
linear slope was subtracted for clarity. The sharp steps in the 
local magnetization in (a) and (c) signify a FO transition. In 
(b), a break in the slope signifies a SO transition. 



of B — H and dB/dH, respectively. A break in the slope 
of the induction is clearly resolved along both ascending 
and descending field sweeps and the derivatives display 
a step structure indicating a SO transition. 

Yet, our main finding is the recovery of the FO ther- 
modynamic phase transition at lower temperatures. This 
is clearly visible in Fig. [TJ: as a reversible discontinuity in 
magnetization B(H) and accordingly, in Fig. [^fc a sharp 
peak in the derivatives dB/dH. This recurrence of the 
FO nature of the transition has been predicted theoreti- 
cally [13] j but has never been observed experimentally. 

These findings can be explained by the following com- 
parison of energy scales. At high temperatures thermal 
fluctuations are dominant enough to weaken the effec- 
tive pinningpotential due to CDs [H, El, E El, El 



: pinning 

18|, EE EE I22J, |2j|, resulting in a FO transition sim- 
ilar to that found in pristine samples. At intermediate 
temperatures, as correlated pinning gains dominance, the 
randomly distributed CDs alter the equilibrium vortex 
matter state, in addition to enhancing irreversible hys- 
teretic behavior. At these temperatures the melting oc- 
curs at intermediate fields B m > B^ and the pristine 
Bragg glass phase is replaced by a porous solid with crys- 
tallites of interstitial vortices imbedded in the pores of a 
rigid amorphous vortex matrix [Til . EH, EH . The size of 
the nanocrystals within each pore is only several times 
their lattice constant. These dilute nanocrystals seem to 
melt into nanodroplets through a SO transition probably 
since the range of correlations is cutoff by the finite size 
of the pores, consistent with recent numerical simulations 
[H, EE EH, [13]. This is similar to the dense CD limit, 
equilibrium measurements at intermediate temperature where correlated pinning dominates and vortex matter 



14J . Accord- 



ingly, B(H) in Fig. [TJl shows a sharp step and the deriva- 
tives dB/dH in Fig. [2^ display 5-like peaks. At a slightly 
lower temperature thermal fluctuations are weaker, the 
effective pinning potential of the CDs gains strength and 
irreversibility is greatly enhanced, thus masking the un- 
derlying phase transition. We overcame this problem by 
applying the shaking method rendering a reversible mag- 
netization 

0, EH El El- Fi gures [TJd and ^Bp thus show 
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undergoes a second order BoG transition. 

At low temperatures ~ 40K the transition occurs at 
high fields. Each pinned vortex is surrounded by tens 
of interstitial vortices and the solid phase is increasingly 
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FIG. 2: (color online) The derivative of the measured induc- 
tion with respect to applied field, dB/dH, of the — 10G 
sample in presence of shaking, measured upon ascending and 
descending field at constant temperatures of (a) 89K — 9IK, 
(b) 76K - 84K and (c) 40K - 42K. Data are shifted verti- 
cally for clarity. A peak in dB/dH in (a) and (c) signifies a 
FO transition, whereas a step in (b) indicates a SO transition. 



dominated by inter vortex interactions [17|, resulting in 
relatively ordered dense vortex crystallites [3]. Conse- 
quently, their melting into nanodroplets is accompanied 
by a discontinuity in the entropy that reflects the sharp 
difference in ordering of the two phases. The FO nature 
of the transition is thus restored, as in pristine samples. 

Note that at intermediate temperatures, where the 
transition is SO, the step in the derivative is positive, 
AdB/dH > (see Fig. [2b). This is consistent with a 
positive step in the specific heat, AC„ > which was 
found in the high-field SO transition of pristine YBCO 
2^, 3(J. It is in contrast, however, to the negative step 
of the derivative, AdB/dT < 0, at the SO glass transi- 
tion reported in pristine BSCCO Further study is 
required in order to explain the origin of these differing 
signs. 

Interestingly, at intermediate temperatures dB/dH 
takes a more intricate shape; rather than a simple upward 
jump, it first decreases, then jumps up and finally slightly 
decreases again. This structure may be explained by the 
association of the step in the derivative at the SO tran- 
sition with the relative ordering of the two phases across 
the transition line. Then, the decrease in the derivative 
might be interpreted as a slight ordering of interstitial 
vortices prior to the major disordering upon melting. A 
similar feature has been observed in the structure factors 
in Monte-Carlo simulations [20]. It is yet unclear, how- 
ever, how these possible structural changes should affect 



300- 



03 



200 



~~ 

_ i 1 1 °». 


/*> CP b 

/ °»/ 

• 

• 

• 

_ • 

• 

B =10G \ cp 

X 

r -i i t - 


X CP . 

• 

• 

• 

B=20G \ CP 


cp cp d 
• 

• 

• - 

B <f 20G \ CP 



30 50 70 90 
T 



30 50 70 90 
[K] 



FIG. 3: (color online) Phase diagrams of irradiated BSCCO 
crystals showing FO (o) and SO (•) transitions. (a) In 
£?0 =5G sample the transition is FO at all temperatures, 
(b) For B4 =10G a FO-SO-FO behavior is found, (c) The 
SO segment expands in B^ =20G sample, (d) A different 
£?0 =20G sample with lower doping that exhibits a FO-SO- 
FO-SO sequence. The values of B$ are shown by dashed lines. 
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the equilibrium magnetization. 

The location and nature of the transition of the four 
samples was measured and mapped. Shaking enabled 
the detection of a melting transition down to 27K. The 
results are summarized in Fig. [3l The FO and SO re- 
gions are marked by open and solid circles respectively. 
In all four irradiated samples, the location of the melt- 
ing line remained similar to that of the pristine parts 
of the samples, with a maximum at ~ 40K. The 5G 
sample exhibits a FO transition at all temperatures. As 
the density of CDs is increased, more of the B — T phase 
space becomes dominated by correlated disorder and con- 
sequently, a larger portion of the transition line becomes 
SO. The 10G and 20G samples both display two CPs 
with a FO-SO-FO sequence. In the 20G sample the two 
CPs are shifted further apart than in the 10G sample and 
the SO transition spreads both to high and low tempera- 
tures. We therefore suggest that this trend persists until 
a full SO BoG transition line is attained in the dense-pin 
limit. 

Note that the FO transition line of irradiated sam- 
ples persists to the left of its maximum in the inverse 
melting region (Fig. These are the first equilibrium 
magnetization measurements in this region in presence 
of CDs. Figure [3ji shows the phase diagram of a sam- 
ple with B$ = 20G and a slightly lower oxygen doping. 
This sample displays even more intricate behavior. Like 
the other B^ = 20G and 10G samples, it also exhibits 
FO-SO-FO behavior. In addition, this sample reveals a 
third CP at extremely low temperatures on the inverse 
melting side of the transition. Below this CP we observe 
a recurrence of the SO nature of the transition. It can 
be attributed to the fact that the melting field decreases 
in the inverse melting region with decreasing T. As a 
result, nanocrystals contain less vortices and correlated 
pinning regains dominance over elasticity. Consequently, 
the transition becomes SO once again, similarly to that 
at intermediate temperatures. 

It is worth pointing out the qualitative difference be- 
tween the behavior reported here and the SO-FO-SO se- 
quence found in pristine samples of the less anisotropic 
YBCO compound @. In irradiated BSCCO all three 
portions of the transition line separate a solid phase from 
a liquid phase. In pristine YBCO, however, the high-field 
SO portion is believed to separate tw o liq uid phases due 
to the existence of a tri-critical point [30( ■ The low field 
SO portion, on the other hand, separates a solid phase 
from a liquid phase and is believed to arise from the in- 
trinsic correlated disorder in YBCO crystals 3ll. |32|. 

In summary, we present thermodynamic evidence for 
new FO-SO-FO behavior of the melting line in BSCCO 
samples irradiated with a low density of CDs. This un- 
usual behavior is due to close competition between four 
different energy scales in the dilute CD limit. As CDs 
are introduced, the melting transition initially alters its 
order from first to second at intermediate temperatures 



where correlated pinning has a dominant effect. At high 
and low temperatures the transition remains FO despite 
the disordering potential of CDs. At low temperatures 
where the melting field is high, the solid phase is domi- 
nated by elasticity rather than correlated pinning due to 
the increasing number of vortices per CD, and the FO 
transition is retained. The FO nature of the transition 
is likewise preserved at high temperatures, where corre- 
lated pinning is weakened by thermal fluctuations. As 
the density of the CDs is increased, the SO segment of 
the transition line expands both to the higher and lower 
temperatures. In one sample a SO transition was also 
found on the low-temperature side of the inverse melt- 
ing portion, resulting in an even more complex FO-SO- 
FO- SO sequence. The observed nucleation and growth 
process of the SO segments along the original FO transi- 
tion line clarifies the general process of transformation of 
phase transitions with increased disorder. In particular, 
it describes the mechanism that leads to transformation 
of a Bragg glass in the presence of point disorder to a 
BoG at high concentrations of correlated disorder, which 
melts through a single SO transition. 
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